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ABSTRACT 

We have mapped the NGC 2023 reflection nebula in [C ii] and CO(ll-lO) with the heterodyne receiver GREAT on SOFIA and 
obtained slightly smaller maps in 13 CO(3-2), CO(3-2), CO(4-3), CO(6-5), and CO(7-6) with APEX in Chile. We use these data 
to probe the morphology, kinematics, and physical conditions of the C n region, which is ionized by FUV radiation from the B2 
star HD 37903. The [C n] emission traces an ellipsoidal shell-like region at a position angle of ~ —50°, and is surrounded by a 
hot molecular shell. In the southeast, where the C II region expands into a dense, clumpy molecular cloud ridge, we see narrow 
and strong line emission from high- J CO lines, which comes from a thin, hot molecular shell surrounding the [C n] emission. The 
[C n] lines are broader and show photo evaporating gas flowing into the C II region. Based on the strength of the [ 13 C n] F=2—1 
line, the [C n] line appears to be somewhat optically thick over most of the nebula with an optical depth of a few. We model the 
physical conditions of the surrounding molecular cloud and the PDR emission using both RADEX and simple PDR models. The 
temperature of the CO emitting PDR shell is ~ 90 - 120 K, with densities of 10 5 - 10 6 cm -3 , as deduced from RADEX modeling. 
Our PDR modeling indicates that the PDR layer where [C n] emission dominates has somewhat lower densities, 10 4 to a few times 
10 5 cm -3 . 

Key words. ISM: Clouds - Submillimeter: ISM - ISM: lines and bands - ISM: individual (NGC 2023) - ISM: molecules - (ISM:) 
photon-dominated region (PDR) 


1. Introduction 


Located at a distance of 350 pc in the L 1630 molecular 
cloud NGC 2023 is one of the most well-studied reflection 
nebulae in the whole sky. This is where the polycyclic aro¬ 
matic hydrocarbon (PAH) m olecules we re first identified in 
the interstellar medium (Sellgren 1984). NGC 2023 has a 
favorable geometry, which makes it an almost ideal labora¬ 
tory for studying the physical conditions in photon domi¬ 
nated regions (PDRs). HD 37903, the B2 V star illuminat¬ 
ing the reflection nebula, is on the near side of the molecu¬ 
lar cloud, while the shape of the reflection nebula indicates 
that the cloud interface is seen largely edge-on in the south¬ 
eastern part of the nebula, where it expands into a dense 
molecular cloud ridge. 

The PDR emission from NGC 2023 has been studied ex¬ 


tensively. Howe et al. 
ter on the Kuiper Air 


(1991) used a Fabry-Perot spectrome- 


Dorne observatory (KAO) to image the 


* The NASA/DLR Stratospheric Observatory for Infrared 
Astronomy (SOFIA) is jointly operated by the Universities 
Space Research Association, Inc. (USRA), under NASA contract 
NAS2-97001, and the Deutsches SOFIA Institut (DSI) under 
DLR contract 50 OK 0901 to the University of Stuttgart. 

** The Atacama Pathfinder Experiment (APEX) is a collabo¬ 
ration between the Max-Planck-Institut fur Radioastronomie, 
the European Southern Observatory, and the Onsala Space 
Observatory. 

*** Herschel is an ESA space observatory with science instru¬ 
ments provided by European-led Principal Investigator consor¬ 
tia and with important participation from NASA. 


nebula in the 158 p m [C ill li ne with a sp atial resolution of 


55", while(Steiman-Cameron et al.l (fT997) observed the far- 

w 


infrared fine structure emission lines of [01] 63 and 145 /im, 
[C n] 158 /rm, and [Si II] a t 35 / im with a cryogenic grating 
spectrometer. Jaffe et al. (1990), using ground based tele¬ 
scopes, studied the PDR region by observing the CO(2-1), 
CO(3-2), CO(7-6), and C 18 0(2-l) transitions and found 
that there was a bright ridge to the southeast of HD 37903, 
where the CO(7-6) emission coincides with [C n] emission 
and fluorescently excited H 2 emission. This ridge was re- 
solv ed into a lumpy f ilamen t in the C91a recombination line 


f 


by Wyrowski et al. (2000), who used the Effelsberg 100 m 


telescope and the VLA with an angular resolution of 11". 
The nebula has been imaged in vibrationally excited H 2 
emission by Gat ley et al. ( 19871) , 

McCa rtney et al.l ( |1999p , and |M 
bright f 


Field et~ah| (199 4} |l9 98|), 
artini, Sellgren & DePoy 
ght fila mentary H 2 emi s sion e specially 


Q1999D , revealing 
to the south and s outhwest. Fleming et al. (2010), Sheffer 
et al. (2011|), and Peeters et al. (|2012|) used the Infrared 
Spectrograph (IRS) on the Spitzer Space Telescope to ob- 


emission m sev- 


serve pure rotational H 2 lines and PAH 
eral regions of NGC 2023 in order to compare the obser ved 
intens i ties w ith predictions from PDR models. Mookerjea] 


et al. (2009) analyzed IRAC and MIPS images from the 
Spitzer Space Telescope as well as SCUBA images from the 
James Clerk Maxwell Telescope (JCMT). They found that 
HD 37903 is the most massive member of a cluster with 20 
- 30 pre-main-sequence (PMS) stars as well as at least two 
young stellar objects (YSOs), MM3 and MM4, embedded 
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in the southern and southeastern part of the reflection neb¬ 
ula. 


2. Observations 


2.1. [C n] and CO(11-10) with SOFIA 


The [C n] emission was map ped wit h _the hetero¬ 
dyne receiver GREAT ^Heyminck et al. 2012) onboard 
the Stratospheric Observatory for Infrared Astronomy 
(SOFIA) during three cycle 1 flights: November 1, 2013, 
February 4 and February 5, 2014. On all flights the GREAT 
LI channel was tuned to CO(ll-lO) in the lower sideband. 
Due to a local oscillator failure the L2 mixer (targeting [C 
ii]) was restored to its basic science configuration for the 
February flights, resulting in somewhat lower beam efficien¬ 
cies due to changes in the optics (see below). The flight on 
November 1 2013 provided excellent observing conditions 
at 13.1 km with system temperatures of ~ 2450 K for the 
LI mixer and ~ 3,000 K for L2. In February the LI system 
temperatures were about the same as in November (2500 - 
2600 K), while the system temperatures for L2 were about 
3,000 K. The beam efficiency, for the November 2013 
flight series was measured to be 0.67 for both channels, 
while it was 0.62 for L2 during the February flights. Table [l] 
lists rest frequencies, half power beamwidths (HPBWs) and 
beam efficiencies for all receiver configurations. All beam 
efficiencies are based on observations of Jupiter. 


The whole nebula was mapped in total-power on-the-fly 
(OTF) mode with 1 s integration time at each dump, and 
7.5" step size. The reference position, at +400", —90" off¬ 
set from HD 37903, has some faint (a few K) [C n] emission 
at ^ 8 kms -1 , which was corrected for in the post pro¬ 
cessing. The map was built up by doing a number of small, 
overlapping maps. The final map size is ^ 6.6' in RA and 
~ 6.8' in Dec, covering most of the optically visible reflec¬ 
tion nebula. Long (5 -10 min) integrations were obtained at 
three positions: (0",0"), (+30",+15"), and (+45",—45"); all 
measured relative to HD 37903 (a(2000.0) =5 /l 41 m 38. s 388, 
J(2000.0) = —02° 15' 32.5"). The spectra had very good 
baseline stability and in most cases nothing higher than a 
2nd order polynomial baseline was removed. The final map 
cubes were created using CLAS^] and were re-gridded to 
6" spatial grid for both [C n] and CO(ll-lO). Maps were 
created with several different velocity resolution by resam¬ 
pling the spectra to either 0.25, 0.5 and 1 kms -1 , in order 
to get a better idea of the morphology and velocity struc¬ 
ture of the emission. The rms noise per resolution element 
is ~ 1.6 K for the C(ll-10) map and 2.5 K for [C n], both 
for a velocity resolution of 0.25 kms -1 . 


1 GREAT is a development by the MPI fiir Radioastronomie 
(Principal Investigator: R. Gusten) and the KOSMA/ 
Universitat zu Koln, in cooperation with the MPI 
fiir Sonnensystemforschung and the DLR Institut fiir 
Planet enforschung. 

2 CLASS is part of the Grenoble Image and Line 
Data Analysis Software (GILDAS), which is provided 
and actively developed by IRAM, and is available at 
http: //www. iram.fr/IRAMFR/GILDAS 


2.2. APEX observations of CO(3-2), 13 CO(3-2), CO(4-3), 
CO (6-5), & CO (7-6). 

NGC2023 was observed on December 1, 2013 using 
the FLASH+ receiver on the 12 m Atacama Pathfinder 
Experiment (APEX) telescope, located at Llano de 


Chajnantor in the Atacama desert of Chile (Giisten et 


al. 


2006). FLASH+ is a dual channel heterodyne SIS re¬ 


ceiver operating simultaneously on orthogonal polarizations 
i n the 345 G Hz a nd in the 460 GHz atmospheric windows 
(Klein et al. 2014). Both bands employ state-of-the-art side- 
band separating SIS mixers. The mixers provide large tun¬ 
ing ranges enabling simultaneous observations of 13 CO and 
12 CO(3-2) in the 345 GHz window and 12 CO(4-3) in the 
460 GHz window. The backen ds are Fast Fourier Transform 
Spectrometers ( Klein et al.|2012 ) with a total bandwidth of 
4 GHz and 76.3 kHz frequency resolution. The observations 
were carried out in good observing conditions resulting in 
system temperatures of 150 - 180 K for 13 CO and 12 CO(3- 
2) and 340 - 380 K for 12 CO(4-3), respectively. A region 
of 5.3' x 6.3' was mapped in on-the fly mode with a spac¬ 
ing of 6" in RA and Dec and all scans were repeated once. 
The maps were created using CLASS and calibrated in T m ^ 
using main beam efficiencies measured on Jupiter, appro¬ 
priate for a somewhat extended source (see Table fTI) . The 
quality of the individual spectra were excellent and only a 
linear baseline was removed from the data. Map cubes were 
created with the same velocity resolution and on the same 
spatial grid as for the SOFIA GREAT observations. For a 
velocity resolution of 0.25 kms -1 , the rms noise in the 13 
CO(3-2), CO (3-2) and CO(4-3) are about 0.4, 0.3, and 
0.6 K, respectively. 

On October 31, 2014 the same area was mapped 
in CO(6-5) and CO(7-6) using CHAMP+ on APEX. 
CHAMP+ is a dual channel 7-pixel hexagonal SIS-mixer 
heterodyne array with a central pixel. It works simul- 
taneonsly in the 450 a nd 350 p m atmospheric windows 
(Kasemann et al. 2006). The weather conditions were ex¬ 
cellent with precipitable water vapor (PWV) in the range 
0.54 to 0.35 mm. The system temperatures for CO(6-5) 
were on the average ^ 1200 K. At the very end of the ob¬ 
serving run the system temperatures went up to ^ 1600 
K. For CO(7-6) the system temperatures were ^ 4000 K 
for most mixers, except for two mixers which were about 
2000 K higher. Because the area covered is quite large, the 
field was divided into 10 overlapping subfields. Observations 
were done in OTF mode, sampling data every 4". The ref¬ 
erence position was at —1385", —274" relative to HD 37903 
and has been verified to be free of CO emission. Each field 
was covered twice, first scanning in RA and then in DEC, 
resulting in maps with fairly uniform depth except at the 
outskirts of the coadded map where the noise level is nat¬ 
urally higher. All spectra are calibrated in T m ^ using main 
beam efficiencies measured on Jupiter (see Table fH. The 
pointing was checked on the nearby a Ori in the CO (6-5) 
line and is good to 2.5". For a velocity resolution of 0.25 
kms -1 , the rms levels are ~ 1.8 K and 3 - 3.3 K for CO(6- 
5) and CO(7-6), respectively. 


3. Results 

3.1. Overall morphology 

The [C n] emission appears to trace an ellipsoidal (egg- 
shaped) region oriented approximately southeast to north- 
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Table 1. Observing setup 


Receiver 

Molecular transition 

Frequency 

[GHz] 

@FWHM 

["] 

Vmb 

GREAT LI 

CO(ll-lO) 

1267.01449 

23.0 

0.67 

GREAT L2 

[C II] 2 F 3/2 -+ 2 P 1/2 

1900.53690 

15.3 

0.67 a 

GREAT L2 

[ 13 C II] F = 2-1 

1900.46610 

15.3 

0.67 a 

FLASH+ 

13 CO(3-2) 

330.587965 

18.5 

0.68 

FLASH+ 

CO(3-2) 

345.795990 

17.7 

0.68 

FLASH+ 

CO (4-3) 

461.040768 

13.3 

0.58 

CHAMP+ 

CO(6—5) 

691.473076 

9.1 

0.49 

CHAMP+ 

CO (7-6) 

806.651806 

7.7 

0.48 


7]mb = 0.62 for the February 2014 flights. 


west (position angle r\j — 50°) (Fig.0. The emission is very 
strong in the southeast where the [C II] region expands 
into a dense molecular cloud ridge, which was mapped by 
Wyrowski et al.l (2000) in HCN(l-O) and HCO + (1-0). In 
the northwest the [C il] emission is more diffuse and ex¬ 
tended. In the southeastern quadrant of the nebula the 
emission is very bright with main beam brightness tem¬ 
peratures between 30 - 75 K. The [C n] emission shows the 
same morphology as the PAH emission (at 8 /i m), PACS 70 
fi m emission, and vibrationally excited H 2 emission, while 
C91 a is much more localized and only seen in the south¬ 
eastern part of the nebula where the PAH emission is very 
strong (Fig.[T|) . Both PAHs and vibrationally excited H 2 
are excited by FUV radiation. The 70 /im emission is dom¬ 
inated by emission from small dust grains, which are also 
heated by photoelectric heating from the FUV radiation. 
Therefore it is clear that the excitation of the [C n] in this 
area is also dominated by the FUV field from HD 37903, 
the B2 star illuminating the nebula. 

Although low- J CO emission is detected over the entire 
mapped region (Fig. [2b , the strongest emission is to the 
south and southeast of the nebula. The 13 CO(3-2) emis¬ 
sion, which traces the column density rather than the tem¬ 
perature of the gas, shows the dense ridge to the southeast. 
It also has a secondary peak coinciding with the northern 
[C n] emission peak at (^ 20", 110", ) indicating that the 
[C n] emission is brighter where the FUV radiation illu¬ 
minates a denser region of the surrounding cloud (Fig.j2|. 
The CO(11-10) emission is only visible in the southeastern 
quadrant of the map, where the C II region is bounded by 
the dense molecular cloud located southeast and south of 
HD 37903. The vibrationally excited H 2 emission imaged 
with high spatial resolution indicates that the interface be¬ 
tween the ionized inner region and the surrounding dense 
molecular cloud is far fr om smooth (Fie ld et al.|1994| 1998 
McCartney et al.| [l999), with lumpy ridges and bright fila- 


mentary structures. The same is true for the PAH emission 


(Fleming et al.|2010 Sheffer et al. 12011 Peeters et ah|2012). 

Our maps of [C nj and high-J CO lines (Fig.[2] & [3| look 
smoother, mostly because of insufficient spatial resolution 
to see such details. The [C ll]-map (Fig.Jl]) looks very sim¬ 
ilar to o ne of the e arly vibrationally excited H 2 2.12 /i m 
maps ([G atley et al.||1987| ), imaged with 20"-resolution. 

Inspection of the [C n] and 13 CO(3-2) channel maps 
(Fig.[3] & J4| show that there is a clear velocity gradient 
from south to north. This velocity gradient is very obvious 
in the three position velocity plots (Figure |5|, which cut 
through HD 37903 at three different position angles: east 
to west, southeast to northwest and north to south. The 


systemic velocities are ~ 10 km s -1 in the south, ^ 10.5 
km s -1 in the southeast and east, while the velocities are 
^12 km s -1 in the north and northwest. The 13 CO(3-2) 
channel maps show a shell-like structure surrounding the 
whole nebula at red-shifted velocities (Fig. [4]), indicative of 
an expansion of the nebula. 


o 

CD 
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“O 
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Fig. 3. Velocity-channel maps for [C n] in gray scale over¬ 
laid with contours at 5 and 10 K, and from there to 70 K in 
steps of 10 K. The ‘+’ shows the position of HD 37903. 


The low-J CO channel maps (Fig. B.l & B.2) show the 
same velocity gradient as seen in [C Ilj and i3 CO, but in 
addition they show blue- and red-shifted emission in the 
southern part of the image, with the strongest high veloc¬ 
ity emi ssion near o r centered on the extreme Class I object 
MM3 (Mookerjea et al. 2009), also known as Sellgren D or 
MIR-63. MM3, w hich is a relatively strong millimeter/sub- 


milli meter source (Wyrowski et al. 2000 Mookerjea et al. 
2009), definitely drives a bipolar molecular outflow, but 


there appear to be several other fainter outflows in the re¬ 
gion, including a blue-shifted outflow lobe from Sellgren 
C (MIR-62) coinciding with the Southern Ridge (SR). All 
outflows are associated with young Class II or Class I ob¬ 
jects. These will be discussed in more detail in a separate 
paper (Sandell et al. 2015, in prep). There is a red-shifted 
velocity feature in the northwestern corner of our CO maps 
without any associated blue-shifted emission. This is most 
likely part of the surrounding molecular cloud, which ap¬ 
pears to be more red-shifted in the north and northwest. 
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nebula. The CO(11-10) emission, however, is only seen in 
the southeastern quadrant of the map, the region with the 
strongest PDR emission. Since the peak brightness temper¬ 
ature of CO(ll-lO ) is < 22 K , which is much less than the 



3.2. A thin, hot molecular shell? 


Jaffe et al. (1990|), who observed the PDR region in 


NGC2023 in CO(2-l), CO(3-2), CO(7-6), C 18 0(2-l), and 
[C n], found that southeast of HD 37903 the CO(7-6) lines 
peak in a bright ridge coinciding with [C II] and fluores¬ 
cent ly excited H 2 emission, and is surrounded by a dense 
shell of warm molecular gas. They concluded that this emis¬ 
sion must originate in a region where the density is ^ 10 5 
cm -3 and where the kinetic temperature is > 85 K. Since 
our observations include completely sampled maps of the 
reflection nebula in both low-J CO lines (CO(3-2) and 
CO(4-3)), as well as in high- J CO lines (CO(6-5), CO(7-6) 
and CO(11-10)) and [C 11 ] with high spatial and spectral 
resolution, we are uniquely positioned to detect where the 
warm, hot, and ionized gas resides in the nebula. 

The star illuminating the reflection nebula, HD 37903, 
is a B2 V star. The low visual extinction toward the star 
suggests that the star lies on the near side of the molecular 


cloud (Howe et al. 1991). The strong PDR emission south¬ 
east of the star looks like a half shell and suggests that here 
the expanding [C 11 ] shell is seen almost edge on. 

In order to understand the geometry of the hot molec¬ 
ular shell, the CO(11-10) emission is first examined. 
This emission only probes hot gas associated with the 
PDR/molecular cloud interface with negligible contribution 
from the cold molecular cloud surrounding the reflection 


100 0 -100 
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Fig. 4. Velocity-channel map for 13 CO(3-2) in grey scale 
overlaid with contours at 2, 5, and 10 K, and from there 
to 50 K in steps of 10 K. The ‘+’ shows the position of 
HD 37903. 


85 K that Jaffe et al.l (1990) inferred from their CO(7-6) 
observations, it must be largely optically thin or sub ther¬ 
mally excited. Therefore it is a good tracer of the hot gas 
in the PDR. The CO(7-6) and CO(6-5) emission is also 
dominated by the hot PDR emission with CO(6-5) having 
T mb ~ 85 K, and CO(7-6) ~ 60 K. The hot PDR emission 
is also seen in CO(4-3) and CO(3-2), but these lines are 
dominated by optically thick emission from the surrounding 
molecular cloud and are strongly affected by self-absorption 
from colder foreground gas. This makes it difficult or impos¬ 
sible for these lines to separate the PDR emission from the 
emission originating in the surrounding molecular cloud. 

The CO(11-10) channel maps (Fig. 6|, especially the 
panel centered at 10.5 kms -1 , indicate the eastern bound¬ 
ary of hot gas to be ~ 90 - 110" offset from the central 
star, HD 37903. This is seen more clearly in the CO(6-5) 
and CO(7 - 6) ch annel maps at velocities 10.5 and 11 kms -1 
(Fig.ro fc |B.4| ), although the boundary is not completely 
smooth. Figure [7] shows three position velocity plots over 
the southeastern quadrant of the map, which cut through 
the molecular cloud ridge in the east, south east and south. 
These cuts show that the CO emission ends at ^ 100" from 
the star, which is assumed to be the radius of the shell (0.19 
pc). The cuts in CO(7-6) and CO(6-5) (Fig. [ 5 ]) show the 
same. In these cuts, however, some faint CO emission is seen 
all the way to the edge of the map, indicating that there is 
some emission from the surrounding cloud even in CO tran¬ 
sitions as high as CO(7-6). At the boundary the CO(11-10) 
emission is centered at 10.5 kms -1 . The CO(ll-lO) lines 
are generally very narrow, ~ 0.5 kms -1 , suggesting that 
the shell is viewed tangentially. The CO(7-6) and CO(6-5) 
lines are broader (>1 kms -1 ), probably because of con¬ 
tribution from the more turbulent surrounding molecular 
cloud. The CO lines get broader or double peaked closer to 
HD 37903 (Fig.[7|. This is to be expected if the line of sight 
goes through the front and the back of an expanding shell or 
a ridge inside the shell. The narrow line widths along the 
boundary of the CO(11-10) emission region suggest that 
the emission resides in a thin shell. Figure [8] shows a sim¬ 
plified picture of the expanding C II region surrounded by 
the hot PDR region. 

Not all the CO(11-10) is emitted from the shell, how¬ 
ever, because there is definitely structure in the CO emis¬ 
sion even at a resolution of 23". This emission appears to 
originate in dense clumps, ridges and filaments inside the [C 
II] region or in the interface between the surrounding dense, 
lumpy molecular cloud and the [C 11 ] region. At the posi¬ 
tion of the CO peak (55", —15") the CO(ll-lO) emission 
appears more like a northsouth ridge with a linear extent 
of 30". The CO lines are broad, 1.5 kms -1 , and are not dou¬ 
ble peaked. This strong CO(11-10) emission region, which 
also stands out in CO(7-6) and CO(6-5), coincides with 
a prominent north south emission ridge, which is bright 
in both PAH emission and vibrationally excited H 2 . The 
northsouth cut goes thr ough t he easter n edge of the south¬ 
ern ridge (SR) (see e.g. S heffer et al.|2 011) at ~ 75" south, 
which stands out as a clear peak both in [C 11 ] and high-J 
CO lines. The emission from the SR is redshifted (Vi 3r ^ 
11 kms -1 ) placing it near the backside of the nebula. 
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Fig. 6. Velocity-channel map for CO(ll-lO) with contours 
at 2.5 to 30 K in steps of 5K. The ‘+’ shows the position 
of HD 37903. 


The cloud in which the reflection nebula is embedded 
is much less dense in the north and northwest, which is 
evident from the much fainter 13 CO emission (Fig.0. Here 
the densities in the PDR shell are too low to excite LJO(11- 
10f| but sufficiently high for CO(6-5) and CO(7-6), which 
are Doth dominated by the hot gas in the PDR (Fig. [2]). In 
the northwestern quadrant one can see double split lines in 
both CO(7-6) and CO(6-5). Here the separation between 
the two velocity features is ^ 2 kms -1 , suggesting that the 
expansion velocity is ~ 1 kms -1 . In the southeastern part 
of the nebula the front and back side of the PDR is barely 
resolved in CO(11-10) and are separated in velocity by only 
~ 1 kms -1 , corresponding to an expansion velocity into 
the dense ridge of ~ 0.5 kms -1 . These expansion velocities 
are consistent with the linear extent of the C II region, 
which is about twice as large in the northwest as it is to 
the southeast. 


3.3. The origin of [C li] and C91a emission 

Although [C n] is detected throughout the mapped region, 
the [C n] emission is dominated by the PDR illuminated by 
HD37903 (Fig.[l]). To the southeast, where the C II region 
expands into the dense lumpy molecular cloud ridge, the [C 


3 For gas at 100 K, the critical density for CO (11-10) is 4 

X 10 3 * 5 cm -3 , while it is 1 0 5 and 7 X 10 4 cm -3 for CO(7-6) and 

CO(6-5), respectively (Yang et al. poio). 


Fig. 7. Position-Velocity diagrams along selected directions 
for CO(11-10). The contour levels for PA=90 and 135° are 
10% to 100% (10% step) of the maximum T m b of 18.4 and 
18.8 K respectively. The contour levels for PA=0° are 20% 
to 100% (20% step) of the maximum T m b of 9K. HD 37903 
is at 0". 


n] emission extends out to ~ 100" from HD 37903. In the 
northwest, where the surrounding molecular cloud is more 
diffuse, the emission from the PDR is fainter and extends 
out to 220" from the star. 

The strongest [C n] emission is detected in the south¬ 
eastern quadrant where it is bounded by the dense molec¬ 
ular ridge (Fig. |9|. The position velocity cut through the 
southeastern ridge (Fig. [5| shows strong, somewhat blue- 
shifted emission inside the molecular shell. The [C n] lines, 
^2-3 kms -1 wide, are broader than the CO(ll-lO) lines 
but roughly coincident in velocity. Inside the shell [C n] also 
shows a fainter, broad, blue-shifted wing extending to ve¬ 
locities as low as 5 kms -1 . The [C n] channel maps (Fig. [ 3 ]) 
show that this blue-shifted emission component is present 
in most of the strong PDR emission region in the south¬ 
eastern quadrant. This blue-shifted emission can easily be 
explained by a photoevaporation flow from the dense PDR, 
suggesting that most of the PDR emission is on the back 
side of the nebula. There is also some fainter (< 15 K) [C 11 ] 
emission to the west protruding out from the hot molecular 
shell (Fig. [ 5 ]). This emission, which is mostly blue-shifted, 
probably originates in the warm surface layers of the molec¬ 
ular cloud, where the densities are lower and the FUV radi¬ 
ation from HD 37903 can directly reach the cloud surface. 
Some of the [C 11 ] emission outside the reflection nebula, ^ 
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5 - 10 K, may be externally illuminated by FUV radiation 
from the large H II region IC434, which is located west of 
NGC 2023. 


OBSERVER 




NW 


Fig. 8. Cartoon showing a simplified picture of hot, slowly 
expanding PDR region surrounding HD 37903. The sur¬ 
rounding cloud is shown in grey with denser gas having 
a darker shading. The hot molecular shell is shown in blue 
and the red shows the [C n] PDR layer. In the southeast, 
where the C II region expands into a dense, lumpy molecu¬ 
lar cloud ridge the expansion is slowed down by the dense 
surrounding cloud, while it has expanded much further to 
the northwest. There are places where the PDR shell has 
broken through the cloud allowing radiation from the star 
to directly illuminate the cloud further out. 


In the northwest the molecular cloud has much lower 
density than the southeastern ridge. The northwestern po¬ 
sition velocity plot (Fig. [5]) shows that there are strong [C 
n]-lines up to ~ 220", which coincide with the PDR ridge. 
Inside this region there are areas where the [C n] and the 
high-J CO lines are double peaked suggesting that they 
come from the front and back side of the expanding [C n] 
shell. The same double split lines are seen even in CO (3- 
2 ), see Fig.jloJ at the offset ~ -50",+53". The 13 CO line is 
dominated by emission from the surrounding cloud and the 
column densities in the hot thin PDR shell are too low to be 
detected in 13 CO. The velocity separation between the two 
velocity peaks (see Table[2]), gives an expansion velocity of 
~ 1 kms -1 , agreeing quite well with the linear extent of the 
C II region in the northwest, ^ 220", which is about twice 



Fig. 9. Integrated intensity map of [C n] (in color) over- 
layed with CO(11-10) and C91a contours in black and 
white respectively. For CO(11-10) the contours drawn are 
at 3, 5, 8, 10, 12, 15, 20, 22, 24, 26 times the value of 
a = 1.3Kkms -1 . For C91a contours are drawn at 5%, 
20%, 30%, 45%, 50% to 100% in steps of 10% of the peak 
intensity of 7.3 10 -3 Jy/beam. The star symbol marks the 
position of HD 37903. 


of its extent to the southeast. Faint [C n] emission is also 
present outside the [C n] shell. The [C n] emission extends 
at least up to 230" from the HD 37903, and is red-shifted 
by several kms -1 . At the edge of our map the velocity of 
the [C n] emission is ~ 12.4 kms -1 . 

The C91a emission mapped by Wyrows ki et al.| (2000) 
is only detected in the southeast (Fig. [9]), where the C II 
region interacts with the dense molecular cloud. The C91a 
emission is also inside the hot molecular shell. The south¬ 
ernmost peak of the C91a ridge, however, is in a region 
where no CO(11—10) emission is seen, nor is there any ev¬ 
idence for a dense gas clump based on our 13 CO data. Even 
though [C n] emission is seen, it is faint and there is no en¬ 
hancement in the emission at the position of the southern 
C91<a peak. 


3.4. Moderately optically thick [C n], evidence from / 13 C 
ll]-emission 

If [ 13 C n] is detected, the ratio of [ 13 C n] to [C n] line inten¬ 
sity, coupled with the 12 C/ 13 C abundance ratio, pr ovides a 
measure of the optical depth of the [C n] emission (Boreiko 
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& Betz 1996). Our GREAT observations cover two of the 
three [ 13 C nj hyperfine lines: the F = 2-1 and F = 1-1, 
offset by +11.2 and +63.2 kms -1 relative to the [C n] ve¬ 
locity, and with relative intensi ties of 0. 625 an d 0. 125 of 
the total [ 13 C il] l ine in tensity (Cooksy, Blake fc Saykally 


1986 Graf et al. 


2012 ] Ossenkopf et al. ' 2013| ). Althougn 


the short integration times did not allow a clear detection 
of the [ 13 C n] F = 2-1 line in individual GREAT spectra, 
the F = 2-1 transition is definitely present even in the av¬ 
erage of all spectra in the map cube. However, in a deep [C 
n] spectrum^] observed toward HD 37903 with HIFI on the 
Herschel Space Observatory , the two strongest hyperfine 
lines of [ 13 C n] , the F = 2-1 at 19 00.466 GHz, and F = 
1-0 at 1900.950 GHz ( Ossenkopf et al.||2013 ), with relative 
intensities of 0.625 and 0.250 respectively, are clearly visi¬ 
ble (Fig.|n). The F = 1 -1 line is too faint to be detected 
in any spectra. Using the fitted peak temperatures for the 
[C n] and [ 13 C n] hyperfine line transitions, the ratio of [C 
II] and the sum of the two [ 13 C II] lines, corrected for the 
relative intensity of the two hyperfine transitions (0.875; 
Ossenkopf et al. 2013), is found to be 32.3. 

To improve the signal to noise of the [ 13 C n] line in the 
GREAT data, all positions in the SE quadrant where the 
[C n] intensity exceeds 35 K were averaged. Figure pT] (right 
panel) shows this averaged spectrum. In this spectrum the 
[ 13 C il] F = 2 -1 transition has a peak temperature of 1.28 
K at a Visr = 21.3 kms -1 , while the peak of the [C il] line 
is 41.7 K at 10.0 kms -1 . The difference in velocity, 11.3 
kms -1 , is in excellent agreement with the expected veloc¬ 
ity separation of 11.2 kms -1 between the [C il] line and the 
F = 2-1 transition of [ 13 C il]. The width of the [ 13 C il] line 
is ~ 1.9 kms -1 , while the [C il] line is about twice as wide, 
3.7 kms -1 . The ratio of the peak temperature between [C 
il] and [ 13 C il], corrected for the relative intensity of the 
F = 2-1 line (0.625) is ~ 20.3. This ratio is very similar 
to what was derived from the spectrum with HIFI towards 
HD37903 (see above). If the [C il] line is optically thin, 
the ratio is expected to be ~ 70, the isotope ratio of 12 C / 
13 C measured for Orion (see e.g., | Ossenkopf et~a l. (2013)). 
The [C il] emission is therefore somewhat optically thick 
with an optical depth of 1-2. All spectra in the northwest¬ 
ern quadrant were also averaged. Even here the [ 13 C il] line 
is present, but only as a 3 -a detection. However, the ra¬ 
tio between [C il] and [ 13 C il] is approximately the same. 
Analyzing other regions in the map give similar results. The 
C II region ionized by HD 37903 is therefore somewhat op¬ 
tically thick in [C il], since the [ 13 C il] is readily detected 
in all regions examined. 

Since the [C il] emission is somewhat optically thick over 
the whole nebula, this can be used to determine the excita¬ 
tion temperate for [C il]. Assuming a beam filling factor of 
one, and an optical depth of 2, the observed brightness tem¬ 
perature at the position of HD 37903, gives a lower limit to 
T ex = 85 K. At the [C il] peak, assuming the same optical 
depth, one gets T ex = 130 K. Overall our observations sug¬ 
gest that [C il] has a similar excitation temperature as the 
CO in the PDR, perhaps slightly warmer, or on the order 
of 90 - 150 K. To improve this estimate a few deep [C il] in¬ 
tegrations are needed providing sufficient signal-to-noise to 
enable a more accurate determination of the optical depth. 


4 This is a 29 minute long load switched observation from the 
Herschel Open Time project OTl_tvelusam_l, retrieved from 
the Herschel science archive. 



Fig. 11. [ 13 C il] spectra calibrated in main beam brightness 
temperature and plotted on a common V/ sr scale of [ 13 C il]. 
On the left is the long integration HIFI spectrum toward 
HD 37903. The F = 2-1 is plotted in black and F = 1-0 in 
red overlaid with Gaussian fits in green. To the right is an 
average of all GREAT [C il] spectra brighter than 35 K in 
the SE quadrant of our map of NGC 2023. The F = 2-1 
line is clearly visible at V/ sr = 10.1 kms -1 , see text. 


3.5. HD 37903 


HD 37903 was found to have faint mid-IR e xcess emission 
fro m a re mnant disk or residual envelope by Mookerje a et| 


al. 


(2009). A “broad” (1.9 kms 1 ) but faint line in CO(11 


TO) is seen toward the star (Table [2). However, the posi¬ 
tion velocity cuts in CO(ll-lO) (Fig. 71, which go through 
the star in three different directions, show no enhancement 
at all at the position of HD 37903. Instead there is a red- 
shifted cloudlet at ^ 11.5 kms -1 , peaking ~ 20" southeast 


and south of the star. Wyrowski et al. 


(2000) found weak 


extended continuum emission at 8 GHz (3.3 cm), with a 
peak southwest of HD 37903 at (—10", —12"), which they 
interpreted as optically thin free-free radiation from gas 
ionized by the B-star. This cloud could be associated with 
the ridge of bright PAH emission which runs from north¬ 
east to southwest just south of the star. There is no clear 
[C il] enhancement at this position, nor in C(11-10) or in 
any of the low- J CO maps, except that they all show blue- 
shifted emission feature south west of the star, which peaks 
roughly at (—10", —18") from the star. A possible scenario 
could be that a clump of gas is being pushed away by ra¬ 
diation pressure from the strong FUV radiation close to 
the star, with the side facing the star completely ionized 
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while the emission further from the star is shielded and 
still molecular. 


3.6. The physical conditions of the cloud surrounding 
NGC2023 

The integrated intensity ratio CO(3-2)/CO(4-3) changes 
very little between 0.9 and 1.1. This suggests that both 
the CO(4-3) and CO(3-2) lines are optically thick at most 
positions in the cloud and hence the brightness tempera¬ 
ture is approximately equal to the gas temperature. Using 
the peak brightness temperatures observed in CO(3-2) we 
estimate the gas temperature of the quiescent cloud out¬ 
side the C II region to be ^ 40 K. If the CO emission is 
clumped the gas temperatures could be somewhat higher. 
At the CO(3-2) emission peak there is a strong contribu¬ 
tion from the hot PDR emission and the observed CO(3-2) 
brightness temperature suggests a temperature of > 75 K. 

In order to further constrain the kinetic temperature, 
density and column density of the gas in the quiescent 
cloud and the hot molecular shell a grid of models b ased 
on the non-LTE ra diative transfer program RADEX (van 
der Tak et al.|2007) have been compared with the observed 
line intensity ratios. Mo del inputs are molec ular data from 
the LAMDA database (Schoier et al. 2005|) and CO colli- 
sional rate coefficients from the work oF Yang et al. (2010). 
RADEX predicts line intensities of a given molecule in a 
chosen spectral range for a given set of parameters: kinetic 
temperature, column density, H 2 density, background tem¬ 
perature and line width. Unlike some of the lighter hydrides 
for which the excitation is dominated by the far- and mid- 
infrared background radiation, the only radiation affecting 
the excitation of CO molecules at (sub)-millimeter wave¬ 
lengths is the cosmic microwave background (CMB). Hence 
a value of 2.73 K is assumed as the background temperature 
for all calculations presented here. The synthetic line ratios 
are calculated assuming line widths of 3kms -1 , which is 
somewhat higher than the observed line widths, but such 
line widths are needed in order to capture the PDR emis¬ 
sion throughout the nebula. 

Figure Jl3] shows contours of the observed CO (3- 
2 )/CO(4-3 )7~C 0(6-5)/CO(4-3), CO(ll-10)/CO(6-5) and 
CO(11-10)/CO(7-6) ratios estimated as a function of i ) 
density (n(H 2 )) and kinetic temperature Ik in, with AT(CO) 
fixed) and ii) N( CO) and Ikm as calculated by RADEX. 
The continuous and dashed contours show the minimum 
and maximum values of the observed ratios respectively. 
While the intensity ratios do not constrain the column den¬ 
sities at all, they uniquely demarcate the kinetic temper¬ 
atures and densities of the clouds which contribute to the 
emission of the CO transitions being considered. As seen in 
the CO(3-2) & CO(4-3) maps primarily trace the 
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Fig. | , _ _ 

quiescent c loud. Based on the C 18 0(2-l) observations by 


Jaf fe et al.| (1990) and assuming that C 18 0 is optically thin, 
the column densities in the surrounding molecular cloud 


densities are adopted and the width of the CO emitting 
layer is assumed to be ^ 1". For a normal CO/H 2 abun¬ 
dance ratio of 10~ 4 N(C O) is ~ 10 17 cm -2 , which agrees 


well with what Jaffe et al. (1990) deduced from their ob¬ 
servations. Since the ratios involving CO(11-10) primarily 
trace the physical conditions of the hot molecular shell, the 
Ikin of the shell Is estimated to be ~ 90 - 120 K for gas 
densities of 10 5 - 10 6 cm -3 based on Fig. 
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3.7. Far-IR line cooling from [C 11 ] and CO 

It is well known that the fine structure lines of [C 11 ] 158 
pm is t he d ominant coo ling line in low-density regions 


(Hollenbach, Takahashi & Tielens 1991) providing ^ 0.3 
% of the total FIR (72 - 196 pm) luminosity in our Galaxy 
(Stacey et al.|1985). Since our observations cover almost the 
whole NGC 2023 reflection nebula, it is possible to deter¬ 
mine the [C 11 ] luminosity for a “typical” reflection nebula. 
Integrating the [C 11 ] emission from 4-18 kms -1 over the 
C II region, an area of 71100 arcsec 2 , yields an integrated 
intensity of 8.86 x 10 -13 W m -2 , corresponding to 3.4L 0 
for a distance of 350 pc. Th is agrees well with the [C 11 ] 
luminosity estimated by Howe et al. (199l]), who estimated 
3.8 L 0 (corrected to a distance of 350 pc) from the area 
they mapped with the KAO. NGC 2023 was observed with 
P ACS as part of the Herschel Gould Belt survey (Andre et 
ah||2010 ), but only at 70 and 160 pm. However, the PACS 
wavelength coverage (60 - 210 pm) (Poglitsch et al. 2010) is 
fairly similar to FIR coverage used by Stacey et al.| Q1985 ), 


are as high as 5 x 10 18 - 10 19 cm -2 . Therefore based on 
modeling with RADEX the Ikm of the quiescent cloud is 
estimated to be 35 - 40 K for densities of 10 5 to 10 6 cm -3 . 

The gas densities in the CO(11-10) emitting PDR shell 
have to be at least 4 x 10 5 cm -3 , the critical density of 
CO(11-10). In the north and northwestern part of the neb¬ 
ula, where CO(7-6) is readily detected but not CO(11-10); 
the densities are somewhat lower, ~ 10 5 cm -3 . To get a 
rough idea of the CO column densities in the PDR these 


allowing an estimate of the flux densities for the same area 
used to determine the [C 11 ] luminosity. Since there is strong 
emission from the surrounding molecular cloud, the inte¬ 
grated emission from the reflection nebula was estimated by 
using a background region to correct for the emission from 
the surrounding cloud. From the PACS images the back¬ 
ground subtracted flux densities were found to be 3200 Jy 
and 3300 Jy for 70 pm and 160 /im, respectively. Assuming 
a dust emissivity index, /3 = 1, the ratio of the 70 to 160 
pm flux densities correspond to a dust temperature of 33.3 
K, which was used to estimate the expected flux density at 
100 /im, 4400 Jy. For such dust temperat ures the color cor¬ 
rection of the PACS bands is negligible (Muller, Okumura 
& K laas|[20TT ), except for the 70 pm band, which requires 
a color correction of ^ 5 %. Therefore the effective band- 
widths from the three PACS photometry band^] can be 
used to compute the FIR luminosity for NGC2023, which 
comes to Lpm = 210 L 0 . For NGC2023 the [C 11 ] luminos¬ 
ity is 1.6% of the FIR luminosity, which sounds quite plau¬ 
sible considering that [C 11 ] is much brighter in NGC 2023 
than in the halos of molecular clouds. 

Since CO(11-10) was readily detected in the dense 
southeastern part of the nebula, one can make a rough es¬ 
timate of how much cooling is provided by CO. The [C 11 ], 
CO, and FIR luminosity are determined by integrating all 
our images over the intense PDR emission in the southeast¬ 
ern and central part of the nebula. This is essentially the 
region w here C O(7- 6) and CO(6-5) emission is very strong, 
(see Fig 
For [C II 


B.4 


fc |B.3| ), and covers an area of 31700 arcsec 2 . 
"aMuminosity of 1.9 L 0 is obtained, i.e., almost 
60% of the [C 11 ] luminosity comes from the dense PDR. 
The emission of CO lines is integrated over the velocity 


5 Herschel PACS ICC Technical Note: The Bandwidth of the 

PACS Photometric System, PICC-CR-TN-044, 2013 
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range 8.9 - 12.6 kms -1 , which should cover all the emission 
from the hot PDR. For CO(ll-lO) the luminosity is 35.0 
x 10 -3 L 0 . For the lower J transitions there is substantial 
contribution from the surrounding cloud. Therefore the in¬ 
tegrated emission was corrected using an adjacent reference 
region. The derived luminosities for the 7-6, 6-5, 4-3, and 
3-2 transitions are: 13.7, 10.8, 6.5 one 3.2 respectively; all 
in units of 10 -3 L©. By estimating the contribution of the 
CO transitions which were not observed, results in a lumi¬ 
nosity of - 0.15 L© for CO transitions up to 11-10. The 
contribution from higher transitions is more difficult to esti¬ 
mate, since there are no observations of any CO transitions 
higher than 11-10. It is therefore assumed that higher J 
CO transitions contribute the same amount as all CO lines 
up to 11-10, which is probably good to a factor of 2. The 
total CO luminosity from the dense PDR region is therefore 
^ 0.3 L©, certainly not more than 0.5 L©. Even though 
cooling from CO is by no means negligible, it is still a factor 
of four lower than the cooling from [C n]. Background sub¬ 
tracted flux densities of the PACS 70 and 160 /im images 
for the southeastern region are 2650 and 2400 Jy at 70 and 
160 /im , respectively. In this region the dust temperature 
is 37 K resulting in a flux density of 3400 Jy at 100 /im. 
The color corrections are now negligible and the FIR lumi¬ 
nosity is 160 L©, i.e., the cooling efficiency from [C n] is 
lower in very dense regions, or only about 1.1% of the FIR 
luminosity. 

These results can be compared to the Orion Bar, an¬ 
other well-studied bright PDR illuminated by a high mass 
star and often considered a pro totypical P D R tem pl ate fo r 
high mass star forming regions. Bernard-Salas et ah] ( 2014| ) 
find that in the dense Orion Bar PDR the [O ij 63 fi m 
lined accounts for 72% of the total line cooling, with [C n] 
contributing less than 18% and CO lines ^ 5%, i.e., the 
contribution from CO is ~ 25% of the [C n] luminosity. 
The Orion Bar is therefore similar to the dense PDR emis¬ 
sion seen in the southeastern part of NGC 2023 in the sense 
that the cooling from CO is still lower than from [C n], but 
not by a huge amount. In contrast cooling from molecular 
lines completely dominate over atomic fine structure lines 
lines in high-mass pro tostars, where cool ing from [C n] is 
completely negligible (Karsk a e^ah| [2014). 


3.8. PDR modeling 

3.8.1. Choice of positions in the PDR 

The close agreement in morphology between the integrated 
[C n] e missi on and fluorescent H 2 and PAH emission (Fig. [lj 
Section 3.1) confirms that dominant part of the [C 11 ] emis¬ 
sion arises from the PDR illuminated by HD 37903. This 
agrees well with the PDR modeling by Kaufman, Wolfire 


& Hollenbach (2006), who showed that at least 85% of the 
[C 11J emission originates in the PDR. We find that the 
high-J CO lines, particularly CO(ll-lO) and CO(7-6), are 
also dominated by emission from the PDR (Section pOl), al¬ 
though there is still some contribution (~10 Kkms” 1 ) from 
the surrounding cloud in CO(7-6) in the dense southeast¬ 
ern ridge. It is difficult to estimate the contribution from 
the surrounding cloud, especially since the densest part of 
the surrounding molecular cloud coincides with the PDR, 
as seen for example from the C 18 0(2-1) map by Jaffe et| 


al. (1990). By integrating over the CO(7-6) emission from 


9 - 12 km s 1 inside the southern ridge and estimating the 


background from a comparison region outside the nebula, 
it is found that about half of the CO(7-6) emission and 
slightly more of the CO(6-5) comes from the surrounding 
molecular cloud. In the northwestern part of the nebula the 
cloud emission is negligible in both CO(6-5) and CO (7- 
6 ). Therefore [C 11 ] and CO(11-10) can be used to explore 
the physical conditions in the southwestern quadrant, while 
results from CO(7-6) and CO(6-5) have to be used with 
caution. For the rest of the nebula these transitions work 
fine, while the lower-J CO transitions are dominated by 
the surrounding cloud. Since the physical conditions vary 
quite of lot, nine positions have been selected for modeling, 
located in different parts of the reflection nebula. These in¬ 
clude: (0", 0"), the position of HD37903; (30", 15"), the 
[C 11 ] emission peak; (53", —13"), the CO(ll-lO) as well 
as C91<a peak; (17", 113"), the northern [C 11 ] peak; (45", 
—45"), a position in the CO(ll-lO) ridge; (—49", 53"), a 
relatively benign position in the northwestern quadrant, 
showing PDR emission from both the front and back side of 
the PDR shell; (23", —60"), centered on the C91a clump 
no. 2, (—13", —79"), the SR position, which is the most 
extensively modeled PDR position in NGC 2023, and (60", 
—60"), peak of the low-J CO emission. The selected po¬ 
sitions are discussed in more detail in Appendix [Aj The 
spectral line profiles of all the selected positions are shown 
in FigflOl and the line parameters derived from Gaussian 
fits (described in the following text) are given in Table |2j 


3.8.2. Determination of line parameters 

Since the spatial resolution of observed images varies from 
24" for CO(11-10) down to 8"2 for CO(7-6), the maps have 
been smoothed by re-gridding all maps in CLASS to the 
same angular resolution and on a common grid. In order 
not to miss out on the considerable structure that has been 
detected in our maps on all angular scales, all maps were 
smoothed to 16'/1, the resolution of the [C 11 ] map, and 
CO(ll-lO) was left at its original resolution (24'/2), since 
CO(11-10) is detected only in the southeastern quadrant 
of the C II region. CO(3-2) and 13 CO(3-2) were also left 
at their native angular resolution, since neither of them are 
good tracers of PDR emission and they are not used in the 
PDR analysis. Smoothing introduces relatively minor ad¬ 
ditional calibration uncertainty, since all spectra were cal¬ 
ibrated in T m fr based on observations of Jupiter, which is 
an extended source with a size of ~ 45" and hence encom¬ 
passes most of the power of the error beam. 

As noted earlier the line emission from the reflection 
nebula NGC 2023 show evidence of emission from the qui¬ 
escent cloud, from the PDR, from the high velocity molec¬ 
ular outflows and also from wings due to photo evaporation 
flows (seen in [C 11 ] spectra). For the analysis of the PDR 
emission most of the emission from the quiescent cloud has 
been excluded by considering only CO(6-5), CO(7-6) and 
CO(11-10) along with the relevant component of the [C 
11 ] emission. Molecular outflows are also excluded except 
for the southern ridge (SR) position at ( — 13",—79"). The 
extraction of the emission primarily due to the PDR com¬ 
ponent was achieved using multi-component Gaussian fits 
to all the observed molecular lines and [C 11 ] for all selected 
positions. 

A stepwise approach was followed to disentangle the 
different velocity components contributing to the observed 
line emission. Since 13 CO(3-2) is a good tracer of the 
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Table 2. Gaussian fits to spectra at selected positions. The CO(7-6), CO(6-5) and CO(4-3) spectra have been convolved 
to the same spectral resolution as [C ii]. CO(ll-lO) upper limits are 3a for a 1 kms -1 wide line. All offsets are measured 
relative to the position of HD 37903. 


Offset 

Tracer 

J T-'mbdv 

Rlsr 

Ah 

J Tmbdv 

Rlsr 

AV 

J Tmijdv 

Vlsr 

Ah 

("") 


(K km s -1 ) 

(km s -1 ) 

(km s -1 ) 

(K km s -1 ) 

(km s -1 ) 

(km s -1 ) 

(K km s -1 ) 

(km s -1 ) 

(km s -1 ) 

(0.0,0.0) 

13 CO(3-2) 

31.6i0.3 

10.25i0.01 

1.59i0.02 








CO(3—2) 

155.4i2.3 

10.23i0.02 

2.52i0.05 








CO (4-3) 

165.8T0.6 

10.22i0.01 

2.55i0.01 








CO(6—5) 

131.3i0.5 

10.25i0.01 

2.23i0.01 








CO (7-6) 

84.7T0.9 

10.24i0.01 

2.OliO.03 








CO(ll-lO) 

4.3i0.2 

10.51i0.05 

1.89i0.13 








[C II] “ 

130.3il.3 

10.36i0.02 

3.45i0.05 

6.3il.01 

9.52i0.04 

1.00i0.09 




(30,15) 

i:j CO(3-2) 

8.1T0.4 

9.09i0.03 

1.06i0.05 

23.9i0.1 

10.45i0.01 

1.74i0.01 





CO(3—2) 

47.7i0.5 

9.42i0.01 

1.30i0.01 

103.7i0.7 

10.77i0.01 

2.32i0.02 





CO (4-3) 

49.1i0.0 

9.45i0.01 

1.26i0.01 

112.0i0.5 

10.78i0.01 

2.23i0.02 





CO(6—5) 

36.6i4.2 

9.49i0.25 

l.lli0.25 

96.2i4.2 

10.66i0.25 

2.00i0.25 





CO (7-6) 

21.8T4.1 

9.62i0.03 

1.07i0.09 

65.6i4.2 

10.56i0.06 

2.04i0.05 





CO(ll-lO) 

0.91T0.4 

9.21i0.10 

0.71i0.29 

6.0i0.4 

10.75i0.06 

1.66i0.15 





[C II ] 6 

159.7T3.9 

9.10i0.02 

2.02i0.04 

51.5i4.1 

11.19i0.06 

1.97i0.12 




(45,-45) 

ia CO(3-2) 

56.0i0.4 

10.09i0.01 

1.34i0.01 








CO(3—2) 

186.6T2.7 

10.39i0.02 

2.44i0.04 








CO (4-3) 

205.9T3.2 

10.41i0.02 

2.41i0.04 








CO(6—5) 

85.8T3.5 

9.85i0.02 

1.36i0.02 

92.3i3.5 

11.00i0.03 

1.55i0.03 





CO (7-6) 

19.4i0.5 

9.70i0.02 

0.97i0.03 

112.9i0.5 

10.54i0.01 

1.93i0.02 





CO(ll-lO) 

10.9il.3 

9.95i0.13 

0.78i0.13 

19.1il.3 

10.69i0.13 

1.76i0.13 





[C II] 

55.0T2.1 

8.33i0.09 

3.41i0.09 

125.2i2.1 

10.42i0.09 

3.27i0.09 




(60,-60) 

13 CO(3-2) 

42.9i2.2 

10.18i0.25 

0.98i0.25 

6.7i2.2 

11.17i0.25 

1.34i0.25 





CO(3-2) 

102.4i4.9 

10.38i0.25 

1.61i0.25 

49.1i4.9 

10.60i0.25 

3.57i0.25 





CO (4-3) 

78.2T2.6 

10.34i0.01 

1.36i0.02 

83.0i2.5 

10.49i0.01 

2.79i0.05 





CO(6—5) 




133.0i0.3 

10.37i0.01 

1.58i0.01 





CO (7-6) 




96.1i0.1 

10.35i0.01 

1.41i0.01 





CO(ll-lO) 




15.7i0.9 

10.43i0.03 

1.03i0.07 





[C II] 

21.0i6.6 

7.51i0.27 

2.OliO.57 

56.4i6.5 

10.02i0.12 

2.26i0.26 




(-13,-79) 

i:j CO(3-2) 

51.2i0.5 

10.03i0.01 

2.OliO.02 








CO(3-2) 

137.2i0.1 

10.07i0.00 

4.09i0.00 

71.4i0.1 

8.51i0.00 

7.44i0.00 





CO (4-3) 

152.2i0.1 

10.06i0.00 

4.09i0.00 

50.5i0.3 

7.75i0.01 

6.65i0.10 





CO(6—5) 

50.9T3.5 

8.49i0.25 

3.72i0.25 

31.1i3.5 

9.42i0.25 

1.27i0.25 

48.9i3.5 

10.97i0.25 

1.47i0.25 


CO (7-6) 

24.6T2.1 

8.19i0.12 

2.88i0.29 

31.9i0.4 

9.55i0.03 

1.37i0.04 

36.Oil.4 

10.85i0.02 

1.30i0.00 c 


CO(ll-lO) 







9.8il.O 

10.71i0.06 

1.30i0.15 


[C ii] 




64.2il.9 

8.85i0.06 

3.32i0.22 

86.Oil.8 

10.61i0.04 

2.36i0.09 

(53,-13) 

i:j CO(3-2) 

47.5T0.5 

10.03i0.01 

1.64i0.02 








CO(3-2) 

46.7T0.7 

10.00i0.00 c 

2.41i0.06 

139.8i0.1 

10.40i0.01 

2.46i0.02 





CO (4-3) 

66.4T0.7 

10.00i0.00 c 

2.40i0.04 

135.7i0.3 

10.45i0.01 

2.26i0.02 





CO(6—5) 

37.Oil.1 

10.00i0.00 c 

2.17i0.03 

140.9il.l 

10.43i0.01 

1.988i0.01 





CO (7-6) 




129.0i0.8 

10.37i0.01 

2.00i0.01 





CO(ll-lO) 




34.7il.8 

10.77i0.04 

1.46i0.09 





[C II] 

50.6T7.8 

7.13i0.15 

2.30i0.28 

126.9i8.2 

10.03i0.09 

2.83i0.18 




(23,-60) 

i:j CO(3-2) 

46.8±0.4 

9.99i0.01 

1.31i0.02 








CO(3—2) 

118.4T0.7 

10.00i0.00 c 

2.94i0.02 

17.1i0.5 

10.35i0.01 

1.10i0.00 c 





CO (4—3) 

115.2T0.7 

10.00i0.00 c 

2.82i0.02 

21.8i0.5 

10.54i0.13 

1.10i0.00 c 





CO(6—5) 

115.8±0.0 

10.23i0.01 

1.92i0.01 








CO(7—6) 

90.3±0.8 

10.21i0.01 

1.75i0.02 








CO(ll-lO) 

15.Oil.2 

10.16i0.04 

1.10i0.10 








[C II] 

112.4i2.7 

10.27i0.03 

2.49i0.09 

50.9il.8 

8.19i0.05 

2.42i0.17 




(-50,53) 

i:j CO(3-2) 

26.7il.O 

10.30i0.04 

2.03i0.10 








CO(3—2) 

35.1i4.8 

10.30i0.00 

3.20i0.00 

18.1i4.8 

9.64i0.25 

1.30i0.00 

37.1i4.8 

11.72i0.25 

1.40i0.25 


CO (4-3) 

19.3il.7 

10.30i0.00 

3.28i0.38 

30.9il.4 

9.78i0.02 

1.52i0.07 

37.Oil.3 

11.71i0.03 

1.35i0.05 


CO(6—5) 

10.2il.l 

10.30i0.00 

3.00i0.00 

24.4i0.7 

9.96i0.02 

1.24i0.01 

21.5i0.1 

11.47i0.02 

1.30i0.00 c 


CO (7-6) 
CO(ll-lO) 

<1.9 



17.5il.5 

9.95i0.04 

1.18i0.10 

13.7il.2 

11.22i0.07 

1.30i0.00 c 


[C ii] 




18.6i4.3 

9.24i0.11 

1.35i0.25 

52.7i5.2 

11.39i0.09 

2.21i0.27 

(17,113) 

i:j CO(3-2) 

30.5i0.4 

ll.10i0.01 

1.81i0.03 








CO(3-2) 

91.2i0.5 

ll.00i0.01 

2.33i0.01 








CO (4-3) 

92.0i0.6 

10.99i0.01 

2.30i0.02 








CO(6—5) 

59.7i0.3 

10.98i0.01 

1.92i0.01 








CO (7-6) 

28.3i0.9 

10.94i0.02 

1.55i0.06 








CO(ll-lO) 

<2.9 










[C ii] 

101.2i3.4 

10.82i0.04 

2.72i0.11 








a Red-shifted wing (3.5 K kms 1 ) ignored in Gaussian fit 
b Both blue- and red-shifted wing ( 11.8 K kms -1 ) ignored in Gaussian fit 


c Not fitted, i.e., kept constant 


quiescent surrounding molecular cloud single-component 
Gaussian fits to the 13 CO(3-2) line was performed to ob¬ 
tain the velocity of the quiescent cloud. Similarly CO (11- 
10 ), if present, and CO(7-6) spectra were used to identify 
the velocity of the hot CO in the PDR by fitting a one- 
or two- (depending on whether the two components arising 


from the front and the back of the PDR are discernible) 
velocity component Gaussians to these spectra. The two 
components often blend in the southeastern quadrant. The 
CO(6-5) spectra were fitted assuming contributions from 
both the quiescent cloud and the PDR, and in most cases 
the best results were obtained by locking the velocity of the 
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quiescent cloud to the value determined from 13 CO(3-2). 
The lower J CO lines, CO(3-2) and CO(4-3) are trickier, 
because especially in the southeastern quadrant they suf¬ 
fer from strong self-absorption by the surrounding cloud. 
Reasonable fits to the low-J CO lines typically required two 
or three components by locking the quiescent cloud compo¬ 
nent. However, due to overlap in the velocity of the emission 
from the quiescent cloud and the PDR, even when a good 
fit has been obtained for the low-J CO lines, the intensities 
may not accurately represent either the cloud emission or 
the PDR emission. Thus both CO(3-2) and CO(4-3) were 
excluded in our comparison of the observations with the 
PDR models. Finally the [C n] line profiles were fitted us¬ 
ing the high-J CO emission as a guideline, although it is 
noted that the [C n] PDR emission is always offset by ^ 
0.2 - 0.5 kms -1 from the hot CO emission and the lines 
are generally broader. Additionally [C n] often show high 
velocity non-gaussian wings from photo-evaporations flows 
which have no molecular counterpart. These were masked 
out out during the fitting and their contribution was de¬ 
termined by taking the difference between the total line 
integral and the sum of the area covered by the Gaussian 
fits. The results of the Gaussian fits are given in Table [2| 


3.8.3. Comparison with plane parallel PDR models 

In order to estimate the physical conditions in the pho¬ 
ton dominated regions (PDRs) contributing to the emis¬ 
sion at the selected positions the observed line intensity 
rat ios are compared with the resu lts o f the m odel for PDRs 


by Kaufman, Wolfire & Hollenbach (2006). The physical 
structure in these models is represented by a semi-infinite 
slab of constant density, which is illuminated by FUV pho¬ 
tons from one side. The model takes into account the major 
heating and cooling processes and incorporates a detailed 
chemical network. Comparing the observed intensities with 
the steady-state solutions of the model, allows for the de¬ 
termination of the gas density of H nuclei, tir, and the 
FUV flux (6 eV < hi/ < 13.6 eV), Go, measured in units of 
the Habing (1968) value for the average solar neighborhood 
FUV flux, 1.6xl0 -3 ergs cm -2 s -1 . As noted already in this 
paper the relative contribution of the cloud and the PDR to 
the total emission depends on the observed tracer and the 
location in the nebula. Since the primary focus is to char¬ 
acterize the PDR in NGC2023, only high-J CO lines like 
(6-5), (7-6), (11-10) and [C n] line intensities correspond¬ 
ing to the velocity components identified to be associated 
with the PDR are being used. In Table [3] the component 
number refers to the first, second or third velocity compo¬ 
nent of Table [2] Appendix [A] discusses in more detail the 
positions and velocity components which were modeled. 

Figure[l4] shows results of comparison of the ob¬ 
served intensity ratios [C n]/CO(6-5), [C ii]/CO(7—6), [C 
il]/CO(ll—10) and CO(ll-10)/CO(7-6) with the predic¬ 
tions of the PDR models. It is assumed that the beam-filling 
factors of the PDR components of the different tracers used 
are identical and hence intensity ratios do not depend on 
the beam filling. Each panel in Fig.[l4| shows the observed 
intensity ratios of different tracers as contours in different 
colors. The strength of the FUV radiation in the region has 
been d erived using the same method as described by |Sheffer| 


4.0xl0 3 r©( 

Parravano, Holenbach & McKee 

2003) is used, 

which is appropriate for a B2 V star of 9 M© ( 

Trundle et al. 


Table 3. Results of comparison of observed intensity ratios 
with PDR model calculations. Positions are given as offsets 
in arcseconds relative to HD 37903. 


Position 

Component # 

n H 

10 4 cm -3 

Gg 

10 3 (Habing Field) 

(0,0) 

1 

5.2 

175 

(30,15) 

1 

1.2-5.2 

24.7 

(30,15) 

2 

8.8 

24.7 

(45,-45) 

1 

2.7-49.8 

7.2 

(45,-45) 

2 

3.0-14.2 

7.2 

(60,-60) 

2 

11.7 

4.1 

(-13,-79) 

3 

1.7-23.1 

4.2 

(53,-13) 

2 

4.5-19.0 

9.2 

(23,-60) 

1 

3.9-13.6 

7.2 

(-50,53) 

2 

5.2 

5.3 

(-50,53) 

3 

2.0 

5.3 

(17,113) 

1 

2.1 

2.2 


a Estimated from stellar parameters as explained in the text. 


2007; Hohle, Neuhauser & Schutz]2010) and the distance to 
NGC 2023 to taken to be D=350pc. In each panel the esti¬ 
mated FUV field strength is marked by dashed horizontal 
lines. Since projected distances to the various positions in 
the cloud are being used, the estimated FUV field will be 
an upper limit and varies between 2.2xlO 3 and 1.7xl0 5 in 
units of Habing field. There is no way to know how close the 
PDR interface is at the position of HD 37903, but it is prob¬ 
ably no closer than 0.025 pc (i.e. 15" radius). Overall, all 
the ratios constrain the densities to values between 10 4 and 


10 4 ' 5 cm 3 corresponding to the range of estimated FUV 


intensities. The ratios involving CO(11-10) tend in general 
to indicate somewhat higher densities of 10 5 cm -3 or more. 
Among the positions which were modeled, the posi- 


tion at (—13",—79") ( 
(SR) that has been mo 

corresponds to the Southern Ridge 

ideled using PDR models by Burton 

et al. (1990b); |Draine & Bertoldi (1996 2000); Steiman- 

Cameron et al. (;1997 

); [Kaufman, Wolfire & Hollenbach 

(2006); Sheffer et al. ( 

2011). Most of these analyses used 


resolved CO spectra. Burton et al. (1990b) found that the 
observational data they had at the time agreed well with 
Go = 10 4 and n = 10 4 cm -3 . Their modeling does repro¬ 
duce our observed [C n] intensities extremely well, but 
over estimate the CO( 7-6) int ensity by about a factor of 


five. Steiman-Cameron et al.| (1997) argued that the [C n] 
emission is likely to be dominated by emission from lower 
density gas surrounding t he P DR based on the_clumpiness 


seen in the H 2 maps. Steiman-Cameron et al. (1997) there¬ 


fore separated the [C 11J emission into two components: a 
dense component (n ~ 10 5 cm -3 ) filling 10% of the beam, 


while a lower density (n ~ 10 3 cm -3 ) fills the rest of the 
beam. A more refined analysis using a merged PDR/H II 
region model predicted that the denser component _cor- 
respond s to Gq = 10 4 and n = 5xl0 4 cm -3 ( |Kaufm an, 


Wolfire & Hollenbach 2006). Our velocity resolved [C II 
observations do not agree with this picture. At the posi¬ 
tion of the SR there are two velocity components in [C 11 ] 
(Table |2j Fig. 10). The red-shifted [C 11 ] component (V/,^ 
= 10.7 Urns -1 ^originates in the dense PDR, see Section [A} 
and most of the emission from the blue-shifted [C II] com¬ 
ponent originates in the PDR on the front side of the C II 
region. Our velocity resolved [C 11 ] observations therefore 
suggest that only < 20% of the [C 11 ] emission is likely to 
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come from low density gas inside the nebula. The densities 
obtained here are between 10 4 and 10 5 cm -3 for the PDR 
emission from the SR depending on which transition of CO 
viz., (6-5) or (11-10) that is used. 


4. Discussion 


NGC2023 has served as a testbed for PDR modeling in 
a num ber of st ud ies. The first PDR modeling was done by 
Black & van Dishoeck (1987) using the near IR observations 
of fluoresc ent vib rationally excite d H 2 emission in the 2 fi m 


window by|Hasegawa et al.|([l987|) at the position 80" south 

thi 


of HD 37903. They found that to reproduce the observed 
H 2 emission line they required a density of 10 4 cm 3 , a 
gas temperature, T^ ~ 85 K, and a radiation field Go = 
300 times that of the Galactic background star light. Since 
then the models as well as observational constraints have 
vastly improved. Models utilizing near-IR fluorescent H 2 
lines at the SR (which includes the position 0”,—80"), find 
gas temperatures of the H 2 emitting region of 700 - 1000 
K, densities of ~ 10 5 cm" 3 , and a radiation field Go ~ 8000 
( Draine fe Bertoldi||1996 


2000 ) 


The advantage of using fluorescent H 2 lines to probe the 
PDR e mission is that there are a large number of H 2 lines , 
see e.g. Burton et al. (1992, 1998); McCartney et al. (1999) 


Takami et al.|( 2000). The ratios of these lines constrain the 
physical conditions, as well as the incident FUV field rather 
well. The line ratios also discriminate between shock or flu¬ 
orescent excitation of the line emission, although there may 
be cases where one can see both fluorescent and shock ex¬ 
cited H 2 . The fluorescent H 2 lines are only emitted from the 
hot gas in the PDR or in shocked gas and have no contam¬ 
ination from low density gas inside the C II region or from 
the surrounding molecular cloud. However, they only probe 
the outer hot layers of a PDR interface, and other tools are 
therefore needed to probe the structure of the PDR. The 
general morpholo gy of a P DR has been described in many 
papers (see e.g., Tielens & Hollenbach 1985| Hollenbach, 


Takahashi & Tielens 1991; Kaufman, Wolfire & Hollenbach 

2006). A PDR illuminated by an early B-star has a hot sur¬ 
face layer of Av ^ 1-2, with atomic H and O, as well as 
C + ; a transition to H 2 and C for Ay > 1 -2, where carbon 


starts to transition into a mixture of C and CO; and for 
Av £ 4 - 6 most of the carbon is in CO, while the oxy¬ 
gen not tied up in CO remains atomic to Ay ~ 10. Here 
the PDR transitions into the cold molecular cloud (10-40 
K), which is probed by CO and other traditional molecular 
gas tracers. The gas in the PDR is heated by photoelec¬ 
tric heating from small dust grains and [C 11 ] is excited by 
collisions with hot electrons and hydrogen. The dominant 
cooling lines are the [C 11 ] 158 /im line and the [O 1 ] line at 
63 /im. 

The hot (1200 - 200 K) intermediate layers are probed 
by [C 11 ], [OI], and [C I]; fluorescent and pure rotational H 2 
lines and high- J CO transitions, while the interior layer of 
a PDR is mostly seen in pure rotational H 2 lines and high- 
J CO transitions. This paper focuses on velocity resolved 
spectroscopy of the [C 11 ] 158 /im line and high-J CO lines 
observed with SOFIA and APEX. The high-J CO lines, es¬ 
pecially CO(11-10) identify the dense PDR regions in the 
southeastern part of the C II region and show that CO re¬ 
sides in a hot (^ 100 K), thin layer surrounding the [C 11 ] 
emission. The close correspondence in velocity (a few tenths 
of kms -1 ) between [C 11 ] and CO(ll-lO) and/or CO(7—6) 


indicate that most of the [C 11 ] emission originates in the 
dense PDR with only a minor part of the emission orig¬ 
inating in low density gas inside the C II region, except 
for some regions in the dense southeastern molecular ridge, 
where there is relatively strong blue-shifted [C 11 ] emission 
without any counterpart in CO. The general agreement be¬ 
tween PDR emission and [C 11 ] emission is qualitatively 
seen in our integrated intensity map of [C 11 ] overlaid on 
other PDR tracers (Fig.Jl), which shows that the [C 11 ] 
emission traces the boundary of the C II region (limb- 
brightened) and the strongest [C 11 ] emission coincide s with 
bright PDR filaments. Our PDR modeling (Section |3.8.3| 
see also Table [ 3 ]), also shows that the observed [C 11 ] emis¬ 
sion can be explained with the [C 11 ] emission coming from 
PDR regions with densities ranging from ~1 x 10 4 cm -3 
- 4 x 10 5 cm -3 , although the radiation field is not very 
well constrained. There are, however, positions in the dense 
southeastern ridge, where the lines are heavily blended and 
the [C 11 ] emission cannot be accurately identified as orig¬ 
inating in the hot PDR shell. The [C 11 ] peak (30", 15") is 
such a case. Here the CO(11-10) emission profile is asym¬ 
metric with the dominant peak at 10.8 kms -1 , while the 
strongest [C 11 ] emission (peak temperature ^ 74 K) is cen¬ 
tered at 9.1 kms -1 , where hardly any CO(11-10) emission 
is seen. At this position the CO(7-6) line profile is symmet¬ 
ric and probably heavily contaminated by the surrounding 
colder molecular cloud. Position (45", —45") is somewhat 
similar. At these positions it is quite possible that more 
than 50% of the [C 11 ] emission comes from low density 
gas inside the C II region. Here it would be useful to have 
observations of CO(8-7) (which can be observed from the 
ground, but only under very good observing conditions), 
or CO(12-ll) and CO(13-12), the latter two can be ob¬ 
served with SOFIA and will have no contamination at all 
from the surrounding cloud. However, an even better way 
to determine what fraction of the ionized gas comes from 
the dense PDR vs. lower density gas inside the C II region 
is to do velocity resolved observations of the [O 1 ] 63 /im 
line, which now can be done with the GREAT H-channel. 
The [O 1 ] 63 /im line has a much higher critical density (4.7 
x 10 5 cm -3 ) compared to [C 11 ] (2.8 x 10 3 cm -3 ) (Meixner 


& Tielens 


1993). If the [C 11 ] emission is dominated by gas 


with densities < 10 4 cm 3 , it will not be seen in [O 1 ]. 


5. Summary and Conclusions 

The NGC2023 reflection nebula was mapped in [C 11 ] 
and CO(11-10) with the heterodyne receiver GREAT on 
SOFIA and slightly smaller maps in 13 CO(3-2), CO(3- 
2), CO(4-3), CO(6-5), and CO(7-6) were obtained with 
APEX in Chile. The [C 11 ] emission was found to trace an 
expanding ellipsoidal shell-like region surrounded by a hot 
molecular shell seen in high-J CO lines. The expansion ve¬ 
locity of the PDR shell measured from the high-J CO lines 
is ~ 0.5 kms in the southeast and ~ 1 kms -1 in the north¬ 
west , which is consistent with the size of the C II region. 
The [C 11 ] emission is dominated by dense (n = 10 4 - 10 5 
cm -3 ) gas in the PDR interface between the C II region 
and the surrounding cold molecular cloud, except perhaps 
in the southeast, where there may be a significant contri¬ 
bution from lower density [C 11 ] emission inside the C II 
region. In the northwest, where the surrounding molecu¬ 
lar cloud is less dense, the densities inside the C II region 
appear to be too low to excite [C 11 ], and all the [C 11 ] 
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emission comes from the PDR shell surrounding the C II 
region. Based on the strength of the [ 13 C n] F=2-l line, 
the [C n] line appears to be somewhat optically thick in 
the PDR shell with an optical depth of a few. The high- J 
CO lines are very narrow (0.5 - 1 kms -1 ) and strong, indi¬ 
cating that the CO emission in the PDR shell comes from 
a thin, hot, molecular shell surrounding the [C II] emission. 
The temperature of the CO emitting PDR shell is ~ 90 - 
120 K, with densities of 10 5 - 10 6 cm -3 , as deduced from 
RADEX modeling. The excitation temperature for [C n] in 
the PDR is similar to the CO in the hot molecular shell, or 
perhaps slightly warner. The [C n] lines are broader 2 - 
3 kms -1 ), indicating that the [C n] emitting region is more 
turbulent. PDR modeling indicates that the densities in the 
[C n] emitting region are also somewhat lower, 10 4 cm -3 
to a few times 10 5 cm -3 which is to be expected, since the 
overpressure in the PDR forces the ionized gas to expand 
into the [C n] region. The [C n] luminosity was found to 
be 1.6% of the total FIR luminosity for the whole reflection 
nebula. In the the dense molecular ridge the fractional lu¬ 
minosity of [C n] is lower, 1.1% while fractional luminosity 
of CO may be as high as 0.3%. 

Our high spatial and spectral imaging of [C n] and high- 
J CO lines like CO(11-10) and CO(7-6) gives us a unique 
insight to the morphology, kinematics and physical condi¬ 
tions of the C II region surrounding NGC 2023. The only 
thing missing is additional high-J CO lines like CO(12-11) 
and CO(13-12), which are unaffected by the colder sur¬ 
rounding molecular cloud. These lines, together with al¬ 
ready existing [C n] and CO(11-10) data would help us 
further constrain the radiation field and physical conditions 
in the PDR shell. Equally useful would be velocity resolved 
spectra of the [O i] 63 fi m line, which can be used to dis¬ 
entangle the fraction of diffuse and dense gas in the C II 
region. All these observations can be done with GREAT on 
SOFIA. 
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Appendix A: Comments on the selected positions agree well with our findings. 


( 0 ", 0 ") The position of HD 37903. Probably a blend of two PDRs, 
although they are too close together in velocity to separated. South 
of HD 37903 the high- J CO lines split into two components, but at 
the position of the star they are close to the cloud velocity. [C II] is 
much broader and it was with two velocity components. [C il] also 
shows a red-shifted wing. This red-shifted emission extends toward 
SW and is roughly coincident with the somewhat blue-shifted ridge 
of CO emission curving out toward HD 37903. 

( 30 ",+ 15 ") This is the position where [C n] peaks, T m ^ ~ 80 K. 
It is close to and pr oba bly associated with the bright fluorescent H 2 
ridge [Field et al.| (1998|) call the Seahorse. Two velocity components 
were fitted to all observed lines, one close to the systemic velocity of 
the molecular cloud, 10.5 kms -1 , and one which is blue-shifted, 9.1 
- 9.4 kms -1 ( Table [ 2 }. 

( — 13 ",— 79 ") This position is centered on the southern ridge (SR), 
one of the most extensively modeled PDRs in NGC2023. It is, 
however, a difficult position to model, because it coincides with the 
blue outflow lobe from Sellgren C (Sandell et al. 2015, in prep.) and 
the [C II] line and high-J CO lines show two velocity components: one 
blue-shifted, and one red-shifted. The SR is definitely associated with 
the red-shifted velocity component (V) sr ~ 10.7 kms -1 ), because 
north-south cuts through the SR show that this velocity component 
coincides with the SR. The SR also stands out in the high-J CO 
channel maps, m ost p ro minen tly in the channel centered on 11 
kms -1 , see Figj6] R3 fc |B.4| It is still seen even in CO(4-3), but 
not as clearly (Fig. B"2). The good agreement in velocity between [C 
II] and CO(11-10) emission suggest that the [C il] emission originates 
from the dense PDR. The blue-shifted emission component has a 
systemic velocity of 10 kms -1 , and almost certainly corresponds to 
a PDR on the front side of the C II region. This PDR is less dense 
than the SR, and not detected in CO(11-10). Although it is clearly 
seen in CO(7-6) and CO(6-5) it is hard to reliably estimate how 
much of the emission is associated with the PDR, since both lines are 
affected by the blue outflow lobe from Sellgren C. The [C 11 ] line is 
more blue-shifted than the CO lines and may include low density gas 
from the C II region. Therefore only the PDR component directly 
associated with the PDR is modeled. It is likely that the intensity 
of the high-J CO transitions is somewhat underestimated since 
our observations are centered on the ridge and the emission only 
fills about half the beam. However, since our PDR modeling only 
use line ratios, the error from ignoring beam filling is relatively minor. 

( 53 ", — 13 ") &: ( 23 ",— 60 ") The position (53", —13") coincides with 
the pe ak of CO (ll 10) e missio n peak as well as with the C91ai clump 
no. 3 (|Wyrowski et al.||2000|), while (53", —13") is approximately 
centered on the C91a clump no. 2. No attempt was made to model 
the C91<a clump no.l, because the CO emission is contaminated by 
the strong red-shifted outflow from Sellgren D. For both positions 
strong [C II] emission is seen at about the same velocity as C91a, 
but at (23", —60") the [C 11 ] line is much broader, 2.5 kms -1 vs. 
0.8 kms -1 for the C91a emission line, probably because more of 
the extended emission associated with the ridge is included and 
because [C 11 ] is not very sensitive to the density of the emitting 
gas. A relatively strong blue-shifted [C 11 ] emission feature is also 
seen at both positions without any counterpart in high-J CO emission. 

( 45 ",- 45 ") & ( 60 ",- 60 ") These two positions probe the strong 
PDR ridge southeast of HD 37903. A “long” integration spectrum in 
[C 11 ] and CO(11-10) was obtained on (45", —45"), which is close to 
the peak position in CO(ll-lO). The position (60", —60") coincides 
with the peak in low-J CO lines. Both positio ns s how emission 
from two blended PDR components in [C II] (Fig.[l0] Table [ 2 }. At 
(45", —45") the red-shifted component dominates the emission in 
both [C 11 ] and CO(11-10). The second component is close to the 
cloud velocity. Here [C il] shows a strong blue-shifted wing making 
it difficult to separate the emission from the PDR from the low 
density [C 11 ] inside the C II region, although PDR modeling was 
done for both components. At (60",—60") our high-J CO spectra 
only show emission from the PDR on the backside of the C II region. 
The blue-shifted [C II] component is almost certainly dominated by 
lower density gas inside the nebula. Here only the red-shifted PDR 
compone nt is m odel ed. The (60",—60") position was also observed 
by |Jaffe et al. (1990) in CO(7-6) with a 34"-beam. Their results 


(— 50 ",+ 53 ") This position is close to the symmetry axis in the 
northwestern part of the nebula. All lines, except 13 CO(3-2) show 
double split lines. This is not due to self-absorption, since there 
is no foreground gas that could absorb the CO(7-6) and CO(6-5) 
emission. At CO(4-3) the emission from the backside of the nebula 
is probably affected by self-absorption near the cloud velocity, 
making the emission somewhat more red-shifted than the higher J 
transitions. The PDR modeling therefore only uses the CO(7-6) and 
CO(6-5) transitions and both PDRs are modeled separately. 

(+ 17 ", + 113 ") This position coincides with the northern 1 3 CO(3- 
2) pea k and close to the crossing H 2 filaments, which |Field et al.| 
(1998} called the “triangle” (located at 16", +129"). Here ah tfie lines 
show only one velocity component at the cloud velocity (Table [ 2 }. 
The triangle was observed by Takami et al.| (T2000) in several H 2 
emission lines in the near-IR with ~ Fabry-Perot spectrometer 
providing spectral resolution, R, of 500 - 2200. This positio n is als o 
close to the position (+33", +105") observed by ([Burton et al. 1 998} , 
who obtained spectra from 1- 2.5 fi m, at 230 - 430 spectral resolution. 


Appendix B: CO channel maps 

Figures |B.l| to |B.3| show the channel maps of CO(3-2), CO( 4-3), 
CO(6-5), and CO( 7-6) . The fields of view of the CO(6-5) (Fig. |B.3} 
and CO(7-6) (Fig. |B.4| ) channel maps were restricted to the south¬ 
eastern and central part of the nebula, in orde r to visu alize t he strong 
CO emission from the PDR better. Figures [B.5| and |B.6| show the 
channel maps for CO(6-5) and CO(7-6) but smoothed to the 1671 
resolution of the [C II] data in order to capture the fainter emission 
in the northwestern part of the map better. 



100 0 -100 

<5R.A. (") 


Fig. B.l. Velocity-channel map for CO(3-2) with contours 
at 2, 5, 10 to 80 K in steps of 10 K. The ‘+’ shows the 
position of HD 37903. 
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Fig. B.2. Velocity-channel map for CO(4-3) with contours 
at 2, 5, 10 to 110 K in steps of 10 K. The ‘+’ shows the 
position of HD 37903. 



100 50 0 -50 

(5R.A. (") 


Fig. B.3. Velocity-channel map for CO(6-5) showing the 
southeastern and central part of the map with full resolu¬ 
tion. The contour levels start at 10 K and go to 100 K in 
steps of 10 K. The ‘+’ shows the position of HD 37903. 
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Fig. B.4. Velocity-channel map for CO(7-6) showing the Fig. B.5. Velocity-channel map for CO(6-5) smoothed to 
southeastern and central part of the map with full resolu- 16'! 1 resolution showing the whole area that was mapped, 
tion. The contour levels are at 15 K, and from 20 K to 100 K The contour levels go from 8 K to 80 K in steps of 8 K. The 
in steps of 10 K. The ‘+’ shows the position of HD 37903. ‘+’ shows the position of HD 37903. 
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Fig. B.6. Velocity-channel map for CO(7-6) smoothed to 
16'/1 resolution showing the whole area that was mapped. 
The contour levels go from 6.3 K to 63 K in steps of 6.3 K. 
The shows the position of HD 37903. 
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100 0 -100 -200 100 0 -100 -200 100 0 -100 -200 100 0 


Fig. 1 . Integrated intensity image of the [C n] line at 158 /im in color. Contours correspond to the tracers mentioned in 
the individual panels. The peak intensities (and contour levels as % of peak) are 8 /im continuum image: 5306MJysr _1 
(2% to 10% (2.5% step), 10 to 50% (5% step), 50 to 100% (20% step)), PACS 70/im 1.32 10 4 MJy/sr (10%, 20%, 30%, 
45%, 50% to 100% (10% step)), C91<a: 7.3 10“ 3 Jy/beam (5%, 20%, 30%, 45%, 5 0% to 100 % in s teps o f 10%), H 2 S(l): 
t§3 10 -2 ° Wc m" 2 /beam (10% to 9 0% in steps of 10%). The C91a image is from Wyrowski et al. ( |2000| and the H 2 S(l) 
image is from Gatley et al.| (1987]). The offsets in all figures are measured in arcseconds relative to HD 37903, which is 
marked by a star symbol. 
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Fig. 2. Integrated intensity image of the [C n] line at 158/im in color. Contours correspond to the tracers mentioned 
in the individual panels. The contour levels are 13 CO(3-2): (3-8), (10-13),(15-18)cr with a = 4.2Kkms _1 . CO(3-2): 
(3-8), 10, 15, 20, 25, 30, 32, 35, 38, 40 and 42 times cr = 5-lKkms- 1 . CO(4-3): (3-8), 10, 12, 15, 18, 21, 24 times 
a = 8.6Kkms -1 . CO(6-5): 3, 5.5, 9.5, 11.5, 13, 15, 20, 25, 30, 35, 40 times cr = 5AKkms~ l . CO(7-6): 4, 6, 11, 15, 18, 
20, 22, 24, 26 times cr = S^Kkras” 1 . CO(ll-lO): 3, 5, 8, 10, 12, 15, 20, 22, 24, 26 times a = l.SKkms- 1 . HD 37903^8 
marked by a star symbol. 
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PA = 90 deg 



5 10 15 

Velocity (kra/s) 



Fig. 5. Position-Velocity diagrams for [C il], CO(4-3), 13 CO(3-2), CO(6-5), and CO(7-6) along directions given by 
position angle, PA of 0, 90 and 135°. The contours are at 10% to 100% (in steps of 10%) of the maximum T m b of each 
plot. For each column we give the maximum T m b from top to bottom. For PA=90° the values are 56.9, 72.8, 24.3, 69.5 
and 54.9 K. For PA=135° the values are 53.5, 80.5, 41.5, 78.2 and 62.3 K. For PA=0° the values are 53.6, 64.3, 29.2, 56.6 
and 43.3 K. HD 37903 is at 0". The position angle is measured counterclockwise from north. 
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Fig. 10. Comparison of [C il], CO(ll-lO), CO(7-6), CO(6-5), and CO(4-3) spectra at selected positions in the NGC 2023 
region. The CO(4—3) spectra are scaled by a factor of 1.5 for clarity. CO(ll-lO) is scaled by a factor of three and filled 
with greyscale. For clarity we omitted CO(ll-lO) for the panel (17,113), because the spectrum is rather noisy and a clear 
non-detection. 
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Fig. 12. Ratios of integrated intensity over the velocity range 9-12 kms 1 . From left to right we show: (a) CO(3- 
2)/CO(4-3), (b) CO(6-5)/CO(4-3), (c) CO(ll-10)/CO(6-5), and (d) CO(ll-10)/CO(7-6). The emission in CO(4-3) 
and CO(3-2) are dominated by the quiescent colder molecular cloud, while the higher J transitions trace the hot molecular 
shell. 
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Fig. 13. Plot of intensity ratios based on synthetic models calculated using RADEX. The contours correspond to intensity 
ratios indicated in the panel as a function of the physical parameters. The ratios are color coded, with the solid line 
representing the low limit and the dashed one the high limit for each ratio as determined from the ratio plots in Fig. [T2) 
The color coding and limits are as follows: black (6-5)/(4-3) 0.9 (solid) - 1.10 (dashed); red (11—10)/(7—6) 0.07 (solid) 
- 0.20 (dashed): green (11—10)/(6—5) 0.05 (solid) - 0.16 (dashed); magenta (3-2)/(4-3) 0.8 (solid) - 1.1 (dashed). The 
upper row corresponds to models for which T^ n and gas density (n(H 2 )) are varied keeping the column density (N( CO)) 
constant to the value indicated in the panel. In the lower row results models in which i n and N( CO) are varied and 
n(H 2 ) are fixed to the values indicated in the panel. 23 
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Fig. 14. Contours of observed intensity ratios plotted on the intensity predictions as a function of hydrogen density tir 
and FUV radiation field (Go) from PDR models, for the selected positions in NGC2023. The horizontal dashed lines 
correspond to FUV estimated using the stellar parameters as described in the text. 
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